Selective cloning of Gaussian states by linear optics 
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We investigate the performances of a selective cloning machine based on linear optical elements and Gaussian 
measurements, which allows to clone at will one of the two incoming input states. This machine is a complete 
generalization of a 1 — » 2 cloning scheme demonstrated by U. L. Andersen et al. [Phys. Rev. Lett. 94, 240503 
(2005)]. The input-output fidelity is studied for generic Gaussian input state and the effect of non-unit quantum 
efficiency is also taken into account. We show that if the states to be cloned are squeezed states with known 
squeezing parameter, then the fidelity can be enhanced using a third suitable squeezed state during the final 
stage of the cloning process. A binary communication protocol based on the selective cloning machne is also 
discussed. 
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I. INTRODUCTION 

Basic laws of quantum mechanics do not allow the gen- 
eration of exactly alike copies of an unknown quantum state 
01 |H [H 01 ■ However, approximate copies can be obtained 
by using devices called quantum cloning machine JH]. The 
first of such devices was studied to deal with qubits and then 
a continuous variable (C V) J^] analog was developed J3, Htl ■ 
Thereafter, CV optimal Gaussian doners of coherent states 
based on two quite different approaches were proposed: the 
one relies on a single phase insensitive parametric amplifier 
H OH], the other, which has been also experimentally real- 
ized, is built around a feed-forward loop flllTl - On the other 
hand, the latter is much simpler than the first one, overcoming 
the difficulty of implementing an efficient phase insensitive 
amplifier operating at the fundamental limit. Since the setup 
of this device is based only on linear components, throughout 
this paper we'll refer to it as linear cloning machine. Ref. lfl2tl 
investigated the performances of the linear cloning machine 
when the input state was a single generic Gaussian state (co- 
herent, squeezed coherent or displaced thermal state) taking 
into account the effect of fluctuation of the input state covari- 
ance matrix, variation in the setups beam splitter ratios and 
losses in the detection scheme. 

The aim of this paper is to show that the protocol used by 
the linear cloning machine to clone a single input Gaussian 
state can be generalized in order achieve the selective cloning 
a state chosen between two inputs. The possibility to select 
one of two states may have useful implementation in binary 
communication systems where the two bits are encoded in 
two quantum states and the goal of the communication is to 
sent the information from one sender to two receivers. We'll 
address this problem in the final part of the paper. 

The paper is structured as follows: in Sec.|ll]we describe the 
selective cloning machine and describe the evolution of the in- 
put states by means of the characteristic function approach. In 
Sec. [Ill] the requirements of selective symmetric cloning are 
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FIG. 1: Selective cloning of Gaussian states by linear optics: the 
two input states Qk, k — 1,2 are mixed at a beam splitter (BS) of 
transmissivity ri. One of the two emerging beams is measured by a 
measurement described by the POVM H v (z) and the outcome z is 
forwarded to a modulator, which imposes a displacement gz on the 
other outgoing beam, g being a suitable amplification factor. Finally, 
the displaced state is mixed with the state £13 at second beam splitter 
of transmissivity ti. The two outputs, q\ and c/2, from the beam 
splitter represents the two clones, which may be made approximately 
equal to either gi or Q2 by changing the gain g from +1 to — 1. 



exploited and the input-output fidelity is studied. Sec. [IV] in- 
vestigate the possibility to enhance the cloning fidelity and in 
Sec.|V]a possible application of the selective cloning machine 
to 1 — > 2 binary communication is proposed. Finally, Sec. [VI] 
closes the paper with some concluding remarks. 



II. THE SELECTIVE LINEAR CLONING MACHINE 

The selective cloning machine based on linear optics and 
Gaussian measurement is schematically depicted in Fig. Q] 
Two input states, denoted by the density operators g^, k = 
1,2, are mixed at a beam splitter (BS) with transmissivity 
T\. A Gaussian measurement with quantum efficiency 77 is 
performed on one of the outgoing beams, the outcome of 
the measurement being the complex number z. According 
to these outcomes, the other beam undergoes a displacement 
by an amount gz, where g is a suitable electronic amplifica- 
tion factor, and, finally, the two output states, denoted by the 
density operators <tl and ^2, are obtained by dividing the dis- 
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placed state using another BS with transmissivity r 2 . When 
Ti = r 2 = 1/2, g = 1, 7? = 1, Q2 = Q 3 = |0)(0| and the 
Gaussian measurement is an ideal double homodyne detec- 
tion the scheme reduces to that of Ref. ifTHl . which was shown 
to be optimal for Gaussian clo ning of coherent states and has 
been investigated in Ref. In the following we 

carry out a thorough description of the selective cloning ma- 
chine using the characteristic function approach. 

The characteristic function Xfc(Afc) = x[f?fc](Ak) associ- 
ated with a Gaussian state Qk of mode k = 1,2,3 (see Fig. |T|i 
reads: 



exp • 



(1) 



where Af. = (x^, yfc) T , (• • • ) T denotes the transposition oper- 
ation, er fc is the covariance matrix, and Xf. = Ti~[gk (x, y) T ] 
is the vector of mean values, x and y being the quadrature op- 
erators x = (a + aJ) and y = —75(0 — af), with a and 

being the field annihilation and creation operator. In turn, 
the initial two-mode state g = Q\ ® g 2 is Gaussian and its 
two-mode characteristic function reads: 



x[g](A) = cxp 
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and A = (Ai, A2). Under the action of the first BS the state 
x[(?](A) preserves its Gaussian form, namely 

xM(A)-xb'](A)=ex P {-lA r < rA-zA T x} , (4) 

where g' = C/bs,i Qi 8> g 2 EA U 1, while its covariance matrix 
and mean values transform as 11 511 : 



er er = S BS1 er Sbs.i 



X X = S BS1 X = (Xi,X 2 ) 
A, B, and C are 2x2 matrices, and 
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(5) 
(6) 

(7) 



is the symplectic transformation associated with the evolution 
operator J7bs,i of the BS with transmission t\. Note that g' 
is an entangled state if the set of states to be cloned consists 
of non-classical states, i.e., states with singular Glauber P- 
function or negative Wigner function lfl6l[l7ll . 

The Gaussian measurement with quantum efficiency rj (see 
Fig.[Tj) is described by the characteristic function 



[n„(^)](A 2 ) = iexp{-§Ajcr M A 2 -iA^X M } 



(8) 



with Xm = \/2 (Rc[z], Im[z]) T and <tm = tm(i). The 
probability of obtaining the outcome z is then given by 



p v (z) = Tr 12 [g'l®Il TI (z)} 
1 



(9) 
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where X [I ® n„(z)](A) = xM(Ai) x[n„(z)](A 2 ), 
x[I](Ai) = 27rd^ 2 )(Ai) and <5( 2 )(C) is the complex Dirac's 
delta function. We also introduced the 2x2 matrix £ = 

B + <T M . 

The conditional state g c of the other outgoing beam, ob- 
tained when the outcome of the measurement is z, i.e., 



TWn„(*)] 

p v (z) 



(12) 



has the following characteristic function (for the sake of clar- 
ity we explicitly write the dependence on Ai and A 2 ) 



x[ec](Ai) = / d 2 A 



x[0'](Ax,A 2 )x[lMz)](-A 2 ) 



= exp 



Pr,(z) 

A - CYT X C T 



(13) 



Ai 



[CS- 1 (l M - X 2 ) + Xx] } . (14) 

Now, the conditional state £> c is displaced by the amount gz 
resulting from the measurement amplified by a factor g. By 
averaging over all possible outcomes of the double-homodyne 
detection, we obtain the following output state: 

g d = [ d 2 z Pv (z)D(gz)g c D^gz), (15) 
Jc 

with -D(C) being the displacement operator. In turn, the char- 
acteristic function reads as follows: 

x[f?d](Ai) = exp{-iAf o-aAx-iAfXd} , (16) 

with er d = A + g 2 T, + g(C + C T ) and X d = Xi + gX 2 . 
The conditional state (fT3T > is then sent to a second BS with 
transmission r 2 (see Fig.Q]), where it is mixed with the Gaus- 
sian state 03, and finally the two clones are generated. Note 
that, in practice, the average over all the possible outcomes z 
in Eq. ( fT5l ) should be performed at this stage, that is after the 
second BS. On the other hand, because of the linearity of the 
integration, the results are identical, but performing the av- 
eraging just before the BS simplifies the calculations. Since 
gd is still Gaussian, the two-mode state gf = g& ® 03 is a 
Gaussian with covariance matrix and mean given by 
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Xr — (Xd, X 3 ) T , 



(17) 
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respectively, which, as in the case of Eqs. © and ©, under 
the action of the BS transform as follows: 



(Tf ~+ a out 
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(18) 
(19) 



where sik and ^ are 2x2 matrices, and Sbs.2 is the sym- 
plectic matrix given by Eq. (0 with t\ replaced by t 2 . Fi- 
nally, the (Gaussian) characteristic function of the clone c^, 
fe = 1, 2, is obtained by integrating over Aft, h ^ k, the two- 
mode characteristic function x[g out ](Ai, A2), where gout = 
^BS,2 Qi ® 03 UL 2 , i.e., 



x[a](A fc ) = — / d 2 Ahx[&ut](Ai,A 2 ) 

= cxp {-f ^4 A fe - iA^ Jf fc I 
The explicit expressions of 3l\ and ^2 are 



jr 2 = % /i — ^ (/iXj + f 2 x 2 ) + ^x 3 . 



with 



/1 = fi(n,T2,g) =Vn +5"v / l - T i 



/2 = h{ri,T2,g) =gy/n- \/l - n 



(20) 
(21) 



(22a) 
(22b) 



(23) 
(24) 



whereas, ^ and ^ can be written in a compact form as fol- 
lows: 

M =r 2 (/f in + / 2 2 c 2 + g 2 c M ) + (1 - r 2 ) <t 3 , (25a) 
£/ 2 =(1 - r 2 ) (AVi + /| cr 2 + g 2 c M ) + t 2 ct 3 . (25b) 



III. SELECTIVE CLONING 

From Eqs. d22l and d25l ) we see that the two outgoing states 
Si and ft are generally different. In this paper we'll consider 
the case in which the clones are equal, therefore, In order to 
make them exactly alike, one have to put t 2 = 1/2 and X3 = 
0: in this case, SE\ = 3£ 2 and six — si 2 . A further inspection 
of Eqs. d22l and ( 1251 ) with t 2 = 1/2, shows that the states 
<,k could be quite different from both the input states, being 
the covariance matrices and the mean values vectors a linear 
combination of the input ones. On the other hand, if f 2 (or 
fx) vanishes, then the Gaussian output states depend only on 
<Ti, Xi (or <t 2 , X 2 ), cr 3 and er M . In the following we'll 
investigate thoroughly this scenario. 

After we have chosen the symmetric outputs setup, i.e., 
t 2 = 1/2 and X 3 = 0, we are interested in removing the de- 
pendence on the state, e.g., g 2 from the output states, namely, 
we want to let f 2 vanish; this is achieved when 



.9 = 9i{n) = \/(l-n)/n, 



(26) 



TABLE I: Selective symmetric cloning (n = T2 = 1/2 and X 3 = 
0): changing the value of the electronic gain from +1 to —1 one can 
choose to clone the state q\ or Q2 respectively. Notice that if g = —1 
a unitary transformation at the output is needed in order to obtain the 
right sign of the amplitude 3£\s, k = 1,2. 
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which gives fx = 

ay ay 



.-1/2 



and leads to 



(2r 1 )~ 1 / 2 X 1 
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C 3 



1-Ti 



<T M 



It is now clear that if the first BS is balanced (n 
obtain 



six 



S£ 2 
si 2 - 



Xx 



crx 



c 3 + cmJ 



(27) 
(28) 
1/2), we 

(29a) 
(29b) 



This is the 1 — * 2 symmetric cloning of the state Qx- This con- 
figuration has been experimentally implemented to optimally 
clone coherent states Jlll[ll]. Notice that ffi(l/2) = 1. 

On the contrary, in order to eliminate the dependence on 
the state Qx one needs (we are assuming again t 2 = 1/2 and 
X 3 = 0): 



g = .92(7-1) = — s/tx/{i - ti) , 

which gives f 2 = —(1 — ti) -1 / 2 and leads to 

I r 1 = ^ 2 = -[2(i-r 1 )]- 1 / 2 x 2 



1 



c 2 + <t 3 



n 



1-Tl 



and if T\ 



1/2 one has 

ay _ ay v 

■X \ — ,J j2 — —Ji. 2 

six = ^2 = o" 2 + - (c 3 + <X M ) 



CM 



(30) 

(31) 
(32) 

(33a) 
(33b) 



As a matter of fact, to obtain the actual symmetric cloning of 
the state g 2 we have to implement an unitary transformation to 
change the phase of the output states as follows: 3&h — ► — SCy^. 
Notice that 32(1/2) = -1. 

The results of this Section are summarized in Table Q] in 
the case of symmetric cloning (n = T2 = 1/2 and X 3 = 0), 
one can select the state to clone simply change the value of 
the gain g from +1 to — 1. 



IV. ENHANCEMENT OF LINEAR CLONING FIDELITY 



The similarity between the input state g k and the clone Qj, 
k,h = 1,2, can be quantified by means of the fidelity 0] 



F(g k ^h) = Tr 



'Qk Qi\/6k 



(34) 
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which, for Gaussian states, reduces to 

F n = F(g k ^ h ) 



^/~Dct[o- k + £if h ] + 6 - Vd 

x exp{-i(X fe - & h ) T {v k + sf h )~\X k - 3? h )} , 

(35) 

where S = 4(Det[er fc ] - i)(Det[^4] - |). Note that for pure 
Gaussian states Det[<Tfc] = \, and in turn 5 = 0. In the case 
of symmetric cloning X k = 3->h, the fidelity d35l l reduces to 



0.7 
0.6 
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0.4 
0.3 
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y/Dct[tT k + £/ h ] + 6 - VS 



(36) 



and the cloning machine is said to be universal because of its 
invariance with respect to displacement of the input states. 

It is a matter of fact that we can now maximize Eq. Pot by 
a suitable choice of the state q% (<ti, <t 2 , and <tm being fixed). 
Without loss of generality we assume that the covariance ma- 
trix associated with q$ has the following diagonal form 



er 3 




with 



2^ + 1 



,2.s 



w 22 



2^ + 1 



-2s 



(37) 



(38) 



i.e., a squeezed thermal state with & mean thermal photons 
and squeezed parameter s. We recall that X% = in order to 
fulfill the symmetric cloning requirements. Now, if 



cm = 



A?! A12 



Ai 



(39) 



FIG. 2: Plot of the fidelity F v ((Tk, &3, <tm) in the case of symmetric 
cloning when q\ is a squeezed state with real squeezing parameter r; 
<T3 is chosen to be the covariance matrix a-g of a vacuum squeezed 
state (solid line) or (To of the vacuum state (dashed line). See text for 
details. We set f] = 1. 
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FIG. 3: Plot of the G(r, rj) given in Eq. J41b as a function of the 
squeezing parameter of the input state for different values of rj, form 
top to bottom: rj = 1.0, 0.75 and 0.5. See the text for details. 



are the explicit forms of the covariance matrices of g k , k = 
1, 2, and of the measurement YL v (z), respectively, then we find 
that the fidelity reaches the maximum for (for the sake of sim- 
plicity we do not report explicitly the dependence of 7SI on 
k, being clear what is the input state g k under consideration) 



_ 1 f4 7ll +A : n 
s=- log ' 



4~~° V4722 + Al s 



sr = , (40) 



i.e., Q3 should be a squeezed vacuum state with covariance 
matrix er 3 = er ¥ = iDiag(e 2s , e~ 2s ). Indeed, such a maxi- 
mization of the fidelity requires the knowledge of 711 and 722- 
The result obtained above generalizes the conclusions given 
in Ref. lfl2ll . The linear cloning machine described in lfl2ll . 
used to perform 1 — > 2 cloning of the state Q\, follows from 
the present scheme choosing Q2 = Q3 = |0)(0|, correspond- 



ing to a 1 = 0-3 



(To 



■I2, and er 



M 



-Jl 



lo, which 



describes the covariance matrix of the double homodyne de- 
tection with quantum efficiency 77. From Eq. (I4D1 we see that 
sending the vacuum into the second BS is the best choice only 
if Qi is a coherent state or a displaced thermal state 1 1 211 (in 
both the cases s = and <r 3 reduces to the vacuum state 
covariance matrix being X3 = 0). On the contrary, when 



cr k is the covariance matrix associated with the squeezed state 



D(a)S(r)\0) 
S{r 



cxp{ hr(a 



a,r) 
2 



where D(a) = exp{c 



it - 



a* a} and 



a 2 )} are the displacement and squeez- 



ing operators, respectively, r being the real squeezing param- 



eter, then 2711 = 27 2 



e and the cloning fidelity is opti- 



mized if £>3 is a squeezed state with squeezing parameter given 
by Eq. d40b . Fig. |2] shows the enhancement of the fidelity in 
the case of squeezed state 1 — > 2 cloning when a suitable 
squeezed vacuum state with squeezing parameter s given in 
Eq. d40b is used instead of the vacuum state as input g 3 (see 
Fig [D- The effect of non-unit quantum efficiency can be seen 
in Fig. [3] where we plot the quantity 



r<( \ Fv( (T l> CT s"! — Fn((Tl, (To, <Tm) 

G ( r > V) = — Tn~ — — L —^ ) ( 41 ) 



Fr){(Ti,<To, (T 



M I 



as a function of r for different values of rj. G(r, rj) expresses 
the relative improvement of cloning fidelity. As it is apparent 
from the plot, one has enhancement of fidelity for any value 
of rj as far as the signals show nonzero squeezing. 
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FIG. 4: 1 — > 2 binary communication assisted by the selective 
cloning machine. The "service provider" provides the communica- 
tion channel by mixing the two states q\ and Q2 at a balanced BS 
and by addressing the outgoing beams to the "sender". The sender 
performs a measurement on one of the beams and displaces the other 
by an amount gz, z being the measurement's outcome and g being 
chosen according to the bit the sender wants to encode. Finally, the 
message is split into two clones by means of a second BS. See text 
for detalis. 
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2 BINARY COMMUNICATION 



In this Section we address an application of the selective 
cloning machine to a 1 — > 2 binary communication protocol. 
The goal is to encode a classical sequence (string) S of two 
classical symbols, e.g., "—1" and into a quantum se- 

quence S' of two quantum states, e.g., Q\ and Q2, eventually 
unknown, and to send it to two receivers, which are interested 
not only in the classical message but also in the quantum states 
encoding it. In this case a cloning machine is necessary to 
generate the copies 1Z± and IZ2 of S' . Let us now assume that 
the sender, which possesses the string S, is not able to gener- 
ate S' himself, so he needs a service provider that provides a 
communication channel based on the states Q\ and Q2. How- 
ever, since the service provider does not know S, the commu- 
nication channel should be independent on on the message the 
sender want to send. In this scenario the selective cloning ma- 
chine (operating in the symmetric cloning regime) presented 
above can be a useful tool. 

The 1 — > 2 communication protocol based on the selective 
cloning machine is sketched in Fig.|4]and can be summarized 
in these steps: 

• the service provider mixes gi and Q2 at the balanced BS 
and addresses the outputs to the sender; 

• the sender performs the double homodyne detection 
onto one of the two beams and displaces the other one 
by an amount gz, z being the outcome of the measure- 
ment and g being chosen according to the entries ±1 of 
S; 

• the displaced beam is divided into the two clones 
Si(z) = = ^ k \z) by means of another balanced 



BS, with k = 1 (k = 2) if g = +1 (g 



-1). 



It is worth noting that the selective cloning machine is now 
operating at a "single shot" regime, namely, each clone is ob- 
tained after a single outcome z of the double homodyne de- 
tection and not after a complete measurement onto a state. In 



turn, each clone actually depends on z. Once the receivers get 
the single clone, they need a strategy to decide if the bit was 
+1, corresponding to Q\, or —1, corresponding to Q2. 

In order to illustrate the protocol, in the following we ad- 
dress the simple case in which 



Qi = Q2 



\a)(a\ , 



(42) 



are coherent states, i.e., = 5I2 and X\ = X2 = 
s/2(a, 0) (for the sake of simplicity we take a as real and 
positive). We recall that the clones of Q2 have the amplitude 
with a 7r phase shift (see Table |I| with respect to input one: 
in this way it is possible to distinguish between ^'(z) and 
^ 2 \z). Note that one has 

U B s,m <8> eat4s,i = l°> <0| ® \V2a)(V2a\ . (43) 

One of the possible strategies to distinguish between ^'(z) 
and <f ' 2 ' (z) is performing a homodyne detection, which is de- 
scribed by the POVM l2lll 



n*(e) 



- (Z/ 2 J } \ 



where a 2 = (1 — e)/(4e), e is the detection quantum effi- 
ciency, and ITj, = \y)(y\, with 



y 2 ' 2 ^ H n (y) 



\y) = -tit E 



n— v 



(45) 



being an eigenstate of the quadrature operator y = -^(a + a^) 

of the measured mode. In equation d45l l H n (y) denotes the 
n-th Hermite polynomials. Finally, the decision is taken ac- 
cording to the following rule: if x > x => k = 1, otherwise 
k = 2, x being a threshold value. On the other hand, cy>(z) 
and ^ 2 \z) are not orthogonal, and then we have to evalu- 
ate the probability to infer the wrong state, namely, the error 
probability defined as follows: 



H e (z) = -[P Z (2\1)+P z (l\2)] 



(46) 



where P z (h\k) is the probability to infer the state ^ h \z) 
when the actual state was ^ k > (z), h ^ k. In writing Eq. (146b 
we assumed that the two states are sent with the same a priori 
probability p = 1/2. The explicit expressions of P z (2|l) and 
P z (l|2) read as follows: 



P 2 (2|l)= dxTik il \z)U x (e) 



—00 
+00 
dx Tr 



^(z)IL x (e) 



(47a) 
(47b) 



It is easy to see that because of the choice of the states Q\ 
and Q2, the probability H c (z) is minimum when x = 0. The 
average error probability is then given by 



H e (a,ri,e) = I d 2 zp n (z) H c {z) , (48) 
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FIG. 5: Plot of the average error probability H c (a, rj, e) as a function 
of the amplitude a and different values of the quantum efficiencies: 
we set e = 1.0 and, from bottom to top, r\ — 1.0, 0.5, and 0.75. 



1. 
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a 

FIG. 6: Plot of the average error probability H c (a, r), e) as a function 
of the amplitude a and different values of the quantum efficiencies: 
we set rj = 0.75 and, from bottom to top, e = 1.0, 0.5, and 0.75. 



where p v (z) is the double homodyne detection probability 
given by Eq. (0. In Figs.[5]and[6]we plot Eq. (|48T i as a func- 



tion of the amplitude a and different values of rj and e. As one 
may expect, in order to reduce the error probability one has to 
encrease the amplitude of the input coherent states. 

It is worth mentioning that when gi and Q2 are non-classical 
state, then t/BS.igi ® 62U^ S 1 is entangled and such 

correlations can be used to reveal the presence of an eaves- 
dropper along the communication line by means a suitable 
non-locality test JHHl. 

VI. CONCLUDING REMARKS 

We have addressed the performances of 1 — > 2 selective 
cloning machine based on linear optics and Gaussian mea- 
surement, which allows to clone one of two incoming input 
states. We have shown that this is achieved simply changing 
the gain of a feed-forward loop. Moreover a third Gaussian 
state can be used in the final stage of the cloning process in 
order to enhance the input-output fidelity. We have found that 
for coherent or thermal states this state reduces to the vac- 
uum state, whereas a vacuum squeezed state depending on 
the squeezing parameter of the inputs and on the measure- 
ment should be considered when the states to be cloned are 
squeezed states. Finally, a protocol for 1 — ► 2 binary com- 
munication involving the selective cloning machine has been 
proposed and the average error probability has been evaluated 
for a particular choice of the involved states. 
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